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l.  INTRODUCTION 

The  radar  reflectivity,  Z  ,  measured  by  conventional 
meteorological  radars  does  not  provide  any  indication  of  the  type 
of  hydrometeor  present  in  the  cloud,  and  is  consequently  unable 
to  differentiate  between  liquid  raindrops  and  the  various  forms 
of  solid  precipitation.  It  can  give  no  direct  indication  of 
whether  damaging  hail  is  present  within  a  storm. 

It  is  customary  to  interpret  the  magnitude  of  the 
reflectivity,  Z,  in  terms  of  a  rainfall  rate,  R.  Various 
empirical  relationships  exist  between  Z  and  R,  but  in  any 
particular  storm  the  rainfall  rate  derived  from  z  alone  may  be  m 
error  by  up  to  50%.  This  error  arises  because  Z,  for  raindrops, 
is  proportional  to  ND®,  where  N  is  the  concentration  of  particles 
of  diameter  D,  summed  over  all  drop  sizes.  The  same  Z  car.  arise 
from  different  drop  size  distributions  which  would  have  very 
different  rainfall  rates.  For  example  the  same  Z  would  be 
observed  for  one  particle  of  size  one  cm  and  10^-2  drizzle  drops 
of  diameter  0.1mm. 

The  differential  reflectivity  measurement,  ZDR,  is  able  to 
remove  many  of  these  ambiguities.  It  provides  an  indication  of 
hydrometeor  type  and  is  able  to  provide  a  measure  of  the  raindrop 
size  distribution  and  so  leads  to  more  accurate  estimates  of 
rainfall  rate.  This  information  is  of  use  in  predicting  the 
attenuation  of  communication  links  by  precipitation,  for 
identifying  potentially  hazardous  regions  of  clouds,  and  appears 
to  have  some  predictive  value  in  identifying  clouds  which  are 
likely  to  grow  explosively  in  the  space  of  a  few  minutes. 

In  essence,  the  technique  involves  comparing  the  radar 
reflectivity  (ZH)  measured  using  horizontally  polarised 
radiation  with  the  reflectivity  (ZV)  sensed  using  vertically 
polarised  radiation.  The  differential  reflectivity,  ZDR,  is 
defined  as: 

ZDR  =  10  log  (ZH/ZV). 

Small  raindrops  are  spherical  and  so  have  equal  values  of  ZH  and 
ZV  and  give  rise  to  zero  values  of  ZDR.  As  the  raindrops  become 
larger  they  become  increasingly  oblate  and  so  ZH  becomes  larger 
than  ZV,  ZDR  (in  dB)  becomes  positive,  and  its  magnitude  is  a 
measure  of  the  mean  size  of  the  raindrops.  The  Z  and  ZDR 
information  can  then  be  used  to  estimate  the  raindrop  size 
distribution,  and  so  provide  a  more  accurate  estimate  of  the 
rainfall  rate  than  is  possible  with  Z  alone. 

The  interpretation  of  the  ZDR  of  ice  is  rather  more 
complicated  because  the  shape  of  ice  particles  is  not  a  unique 
function  of  their  size.  If  the  ice  particles  tumble  as  they  fall, 
as  is  usually  the  case  within  convective  clouds,  then  although 
for  an  individual  particle  ZH  and  ZV  may  be  unequal,  on  the 
average  the  two  reflectivities  are  the  same,  the  value  of  ZDR  in 
areas  colder  than  freezing  is  usual  lly  close  to  zero.  In 
stratiform  clouds  the  situation  is  rather  different.  Values  of  Z 
are  lower  than  in  convective  clouds,  but  some  regions  of  positive 
ZDR  are  observed.  This  positive  Z  DR  is  believed  to  be  due  to 
small  single  crystals  of  ice,  such  as  needles  or  plates,  which 
have  their  horizontal  diameter  greater  than  the  vertical  and 
which  fall  without  tumbling. 
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During  the  first  year  of  this  researach  project  the  principal 
task  has  been  the  analysis  of  2  and  ZDR  measurements  made  in  both 
convective  and  stratiform  clouds.  The  appendix  coi  tains  three 
papers  which  describe  the  work  in  detail,  one  is  already 

published,  a  second  is  to  be  presented  at  the  forthcoming  Radar 
Meteroolgy  conference,  and  the  third  has  been  submitted  for 
publication.  The  subject  of  the  first  paper  is  the  use  of  Z  and 
ZDR  for  identification  of  hail  and  is  discussed  in  Section  2. 
The  precise  shape  of  raindrops  of  various  sizes  is  of  crucial 
importance  in  interpreting  the  ZDR  of  rain  in  terms  of  a  mean 
raindrop  size;  an  exhaustive  statistical  study  of  the  variation 
of  Z  with  ZDR  suggests  that  raindrops  are  rather  more  oblate 
than  the  previously  accepted  shapes.  This  work  is  described  in 
Section  3  and  appendix  B.  In  Section  4  of  this  report  (and  in 
greater  detail  in  the  paper  in  appendix  C)  we  analyse  the 
evolution  of  convective  clouds.  We  conclude  that  the 
precipitation  in  vigorously  growing  convective  clouds  appears  to 
grow  via  the  coalescence  of  supercooled  droplets  and  not  via  the 
ice  phase.  Such  clouds  contain  unexpectedly  large  supercooled 
raindrops  which  could  be  hazardous  to  aviation,  and  would  also 
cause  anomalously  high  attenuation  of  em  radiation. 

2  HAIL  DETECTION 

Hail  is  one  of  the  most  dangeroues  aspects  of  convective 
storms.  In  the  past  several  radar  techniques  have  been  suggested 
to  identify  hail,  but  all  have  been  found  to  be  unsatisfactory. 
The  dual  polarisation  radar  technique  appears  reliably  to  be 
able  to  identify  hail  within  clouds.  The  brief  note  to  Nature  in 
appendix  A  describes  the  technique.  In  essence  hail  only  occurs 
in  fairly  vigorous  showers  having  reflectivities  above  50dBZ. 
Near  to  the  ground  the  precipitation  is  usually  in  the  form  of 
heavy  rain  which  would  have  a  fairly  large  mean  drop  size;  these 
large  oblate  drops  would  result  in  an  observed  value  of  ZDR 
always  exceeding  ldB.  However  hailstones  tumble  as  they  fall  and 
and  thus  result  in  a  zero  ZDR  value.  Both  heavy  rain  and  hail 
will  be  associated  with  high  values  of  Z,  but  any  hail  should 
give  rise  to  zero  ZDR  and  be  clearly  distinguishable  from  the 
surrounding  heavy  rain  which  has  positive  ZDR.  This  is 
illustrated  in  figure  1,  which  is  a  vertical  section  in  Z  and  ZDR 
through  a  vigrous  convective  cloud.  Below  2km  altitude  the 
values  of  Z  are  generally  positive  when  Z  is  above  30dBZ,  as 
would  be  expected  for  any  reasonable  drop  size  distribution. 
In  contrast  to  this  is  a  region  at  87km  range  where  zero  values 
of  ZDR  extend  right  down  to  the  ground.  If  this  was  due  to  small 
spherical  raindrops  then  the  computed  rainfall  rates  would  be 
several  meters  an  hour,  and  the  only  plausible  explanation  .  s 
that  the  particles  are  tumbling  hailstones,  which  are  so  large 
that  they  reach  the  ground  without  melting  appreciably. 

2.  PROPOSED  NEW  RAINDROP  SHAPES 

The  shape  of  raindrops  is  critical  in  the  interpretation  of 
the  ZDR  data.  Gans  scattering  theory  is  used  to  calculate  the 
values  of  ZH  and  ZV  expected  for  oblate  spheroids  having  various 
axial  ratios.  The  Chilbolton  radar  is  able  to  measure  ZDR  to 
O.ldB  which  is  equivalent  to  a  change  in  axial  ratio  of  about 
0.01.  This  change  is  less  than  the  error  in  the  best  laboratory 
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measurements  of  axial  ratios  by  Pruppacher  and  Pitter.  Earlier 
work  at  Chilbolton  had  suggested  a  small  change  in  the  proposed 
drop  shapes  of  less  than  the  errors  of  the  laboratory 
measurements,  so  that  the  Z  and  ZDR  observed  with  the  radar  would 
agree  with  the  drop  sizes  measured  directly  at  the  ground  7 u s t 
200m  below  the  radar  beam  and  remove  a  small  sytematic  error.. 

However  there  is  considerable  uncertainty  as  to  the  precise 
shapes  of  drops  having  an  equivalent  diameter  above  4mm.  Direct 
verification  is  difficult.  In  the  laboratory  the  precise  shape 
and  degree  of  oscillation  of  such  large  drops  will  vary  wth  the 
degree  of  turbulence  in  the  wind  tunnel.  In  the  field  the  number 
of  occasions  when  such  heavy  rain  falls  directly  over  the  ground 
based  measuring  system  is  small. 

We  have  adopted  a  statistical  approach.  It  is  known  that  on 
the  average,  the  concentration  of  raindrops  (N)  of  size  D 
follows  the  Mar shall-Pa llmer  size  distribution: 

N  (  D )  =  No  exp  (-3.67  Do/D) 

where  No  and  Do  are  constants.  By  summing  the  contributions  of 
the  various  differently  sized  drops  to  ZH  and  Z V  we  can  compute 
how  we  expect  ZDR  to  increase  as  the  rainfall  rate  increases  and 
Z  becomes  larger.  This  is  shown  in  Figure  2  where  the  dotted 
curve  marked  P  is  for  the  predicted  average  variation  of  Z  with 
ZDR  using  the  widely  accepted  raindrop  shapes  of  Pruppacher  ad 
Pitter.  According  to  these  curves  any  value  of  ZDR  above  3.5dB 
cannot  be  interpreted  in  terms  of  raindrops  and  we  must  invoke 
some  very  distorted  half  melted  hail  pellets.  A  very  limited 
series  of  measurements  made  by  Cooper  et  al  using  2D  probes  to 
measure  the  shape  of  larger  drops  suggested  a  deviation  from  this 
form,  and  the  dotted  curve  C  has  been  computed  by  extrapolating 
their  limited  measurements.  The  solid  curve  is  the  average  of 
90000  data  points  observed  by  the  radar  on  eight  different  days 
during  the  summer  of  1983.  This  seems  to  demonstrate  that  indeed 
the  larger  raindrops  are  much  more  oblate  than  suggested  by  the 
laboratory  experiments. 

There  are  several  implications  for  the  radar  work.  First, 
values  of  ZDR  of  up  to  6dB  can  be  explained  in  terms  of  raindrops 
instead  of  needing  to  invoke  extraordinarily  distorted  melting 
particles.  Second,  the  rainfall  rates  and  drop  concentrations 
in  heavy  thunderstorm  type  rain  are  quite  dependent  upon  the 
precise  shape  of  the  larger  drops,  for  values  of  Z  near  60dBZ 
the  two  models  predict  raindrop  oncen t r a t i ons  differing  by  an 
order  of  magnitude.  Finally  these  more  oblate  drop  shapes  will 
lead  to  different  attenuation  on  microwave  links  than  that 
predicted  by  the  previous  models.  A  nore  detailed  description  is 
given  in  appendix  B. 

4.  STUDIES  OF  CONVECTIVE  PRECIPITATION  DEVELOPMENT. 

An  analysis  of  the  rapid  development  of  convective  clouds 
observed  in  1983  has  been  completed  and  the  rather  lengthy  paper 
describing  this  work  is  included  as  appendix  C.  Here  we  summarise 
the  conclusions  from  this  work: 

a)  Isolated  first  echoes  having  low  Z  but  accompanied  by 
large  positve  differential  reflectivity  (ZDR)  are  likely  to 
intensify  very  rapidly. 
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b)  These  anomalous  low  2  high  2  DR  firs':  echoes  indicate 
the  presence  of  very  low  concentrations  of  supercooled  raindrops 
having  diameters  of  several  millimeters,  which  have  grown  in  the 
absence  of  ice.  Present  theories  of  precipitation  development  do 
not  predict  the  appearance  of  such  large  supercooled  raindrops. 

c)  Such  low  concentrations  of  large  drops  could  be  very 
efficient  hail  embryos  growing  to  large  hailstones  after  freezing 
because  of  the  lack  of  competition  for  the  cloud  wter.  This  may 
relax  the  condition  for  embryo  sorting. 

d)  Columns  of  positive  ZDR  extending  above  the  zero  degree 
isotherm  occur  quite  frequently  during  vigorous  corvection.  They 
coincide  with  periods  when  the  echo  is  moving  upwards  at  about 
lOm/s  and  tend  to  precede  high  values  of  2. 

5  PROGRAM  FOR  THE  SOMMER  AND  FALL  OF  1986 

The  radar  is  being  modified  to  measure  additional 
polarisation  parameters.  At  present  (June  1986)  it  is  envisaged 
that  it  will  be  operational  by  the  end  of  July  1986.  One  extra 
parameter  is  the  linear  depolarisation  ratio  ( LDR )  defined  as 
LDR  =  10  log  (ZHV/ZV) 

where  ZHV  is  the  power  received  in  the  horizontal  for 
transmission  in  the  vertical.  In  essence  LDR  provides  a  measure 
of  the  fall  mode  of  the  precipitation  particles.  Raindrops  fall 
with  their  figure  axis  aligned  in  the  vertical  and  so  should  have 
a  very  low  value  of  LDR,  but  ice  particles  which  tumble  will 
have  a  finite  orthogonal  return.  The  second  parameter  is  the 
pulse  to  pulse  fluctuation  in  ZDR,  which  we  call  SDR,  which 
should  give  an  indication  of  the  width  of  the  size  distribution 
of  the  particles.  for  monodispersed  particles  ZDR  does  not 
fluctuate  and  SDR  should  be  zero,  finite  values  of  SDR  suggest  a 
wide  distribution  of  particle  sizes  are  present. 

We  initially  propose  the  following  tasks: 

a)  Further  work  on  tracking  hail  cells  using  both  the  ZDR  and 
the  proposed  LDR  signature. 

b)  Confirmation  in  widespread  rain  that  the  LDR  of  rain  is 
indeed  very  low  and  can  be  used  as  an  additional  means  of 
distinguishing  rain  from  ice. 

c)  Supercooled  raindrops  in  vigorous  convection. 
Confirmation  using  LDR  that  these  particles  are  indeed 
raindrops . 

d)  LDR  measurement  of  the  ice  in  convective  clouds,  ZDR  is 
zero  indicating  tumbling  particles  which  should  therefore  have  a 
measurable  LDR  return. 
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AlixLude  (Km)  Altitude  (Km) 


A)  THE  RADAR  REFLECTIVITY,  Z,  in  dBZ  above  ltnm^m  3 


A  vertical  scan  obtained  on  6  July  1983  at  1336  UT  through  an  intense  convective 

storm.  In  A)  contours  of  Z  are  drawn  every  lOdBZ  with  an  additional  dotted  5dBZ 

contour.  In  B)  the  contours  of  differential  reflectivity,  Z  ,  are  plotted  at  ldB 

DR 

intervals  starting  at  ldB,  together  with  the  dotted  5dBZ  contour  from  A)  which  serves 
to  delineate  the  extent  of  the  echo.  The  positive  values  of  Z  below  2km  altitude 
are  due  to  raindroos  falling  with  their  horizontal  axis  larger  than  the  vertical.  In 
contrast  to  this  at  87km  range  is  a  narrow  region  less  than  1km  wide  where  zero  7 
extends  down  to  the  ground.  It  is  believed  that  this  is  due  to  tumbling  hail  falling 
to  the  ground. 


FIGURE 


For  explanation  see  Text,  Section  3. 
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Hail  causes  considerable  damage  to  property  tod  agriculture  and 
a  simple,  unambiguous  method  of  determining  which  convective 
storms  contain  hail  would  be  of  great  use  for  short-term  weather 
forecasting,  and  would  also  find  application  in  the  study  of  the 
mechanism  of  hailstone  growth.  Previous  radar  methods  of  bail 
detection1-’  have  not  proved  reliable,  but  here  we  report  results 
of  a  new  technique  based  on  the  measurement  of  differential 
reflectivity  using  dual  polarization  radar  and  suggest  an  additional 
criterion  for  bail  recognition.  Hail  is  known  to  occur  in  vigorous 
convective  storms  in  which  tbe  effective  radar  reflectivity,  Z, 
measured  by  conventional  radar  is  high.  If  the  hail  reaches  tbe 
ground  it  is  usually  restricted  to  long  narrow  swathes4,  often  about 
I  km  wide  but  up  to  50  km  long,  accompanied  by  surrounding 
precipitation  in  the  form  of  heavy  rain.  The  differential  reflectivity 
technique  involves  comparing  the  radar  reflectivity  factor  for 
horizontally  (ZH)  and  vertically  (Zv)  polarized  radiation;  the 
comparison  provides  a  measure  of  the  mean  shape  of  the  precipita¬ 
tion  particles.  Close  to  the  ground  the  values  of  Z  are  high  for 
both  bail  and  heavy  rain,  but  because  the  large  raindrops  are 
oblate  and  fall  with  their  horizontal  diameter  greater  than  the 
vertical,  the  value  of  ZH  in  the  rain  exceeds  that  of  Zv-  The 


smaller  regions  where  hail  is  present  can  be  dearly  distinguished 
from  tbe  neighbouring  heavy  rain  because,  although  individual 
hailstones  may  not  be  spherical,  they  tumble  as  they  fall  and  so 
result  in  equal  values  of  ZH  and  Zv. 

For  meteorological  radars  the  power  scattered  back  by  com¬ 
mon  precipitation  particles  is  proportional  to  the  sixth  power 
of  the  particle  size  and,  although  the  particles  may  be  frozen 
or  liquid  and  may  be  noo-spherical,  it  is  usual  to  express  the 
return  as  the  equivalent  radar  reflectivity  factor,  Z,  defined  as 
the  summation  per  unit  volume  of  the  diameter,  raised  to  the 
sixth  power,  of  spherical  water  drops  which  would  scatter  the 
same  power  as  that  actually  measured.  Z  is  usually  quoted  in 
decibels  above  1  mm6  m-3.  A  full  review  of  meteorological  appli¬ 
cations  of  radar  is  provided  by  Browning5 

The  new  method  of  hail  detection  relies  on  a  measurement 
of  the  differential  reflectivity,  ZDR,  defined  in  decibels  as 

^DR  =  1®  logio  (ZH/Zy) 

The  radar  used  in  this  study  is  fully  described  by  Cherry  and 
Goddard6  Situated  at  Chilbolton  in  southern  England,  it  has 
an  antenna  of  25  m  diameter,  giving  a  0.25°  beamwidth  at  the 
10  cm  operating  wavelength.  Pulses  with  horizontal  and  vertical 
polarizations  are  transmitted  (and  received)  alternately  with  a 
pulse  repetition  frequency  (PRF)  of  610  Hz.  The  precipitation 
forms  an  incoherent  target  with  the  values  of  Z  fluctuating  as 
the  hydrometeors  reshuffle  in  space,  but  this  choice  of  PRF 
ensures  adequate  correlation  of  Z  between  sequential  samples 
so  that  ZDR  may  be  estimated  accurately.  Each  range-gated  data 
sample  comprises  a  spatial  average  over  300m  (four  pulse 
volumes),  which  is  time  averaged  over  64  transmuted  pulse 
pairs.  The  worst  case  errors  in  measuring  Z  and  ZDr  are  esti¬ 
mated  to  be  0.7  dBZ  and  0.2  dB  respectively 

A  venical  section  through  a  vigorous  convective  storm  is 


Fig.  1  Contour  plots  of  reflectivity  Z  (a) 
and  differential  reflectivity  ZOR  <6i 
obtained  from  a  vertical  scan  through  an 
intense  convective  storm  at  1 3 :  36  ut,  6  July 
1983.  a.  Contours  of  Z  are  drawn  every 
lOdBZ  above  I  mm®  m  \  with  an  addi-  _ 
tional  dotted  5-dBZ  contour  b.  Contours  E 
of  ZD„  are  plotted  at  1-dB  intervals  star  ng  £ 
at  1  dB  The  dotted  5-dBZ  contour  from  a  ® 
is  included,  to  delineate  the  extent  of  the  ^ 
echo  The  positive  values  of  ZDR  below  2  km  “ 
altitude  are  due  to  raindrops  fallirg  with  ^ 
their  horizontal  axis  larger  than  the  verncal. 

In  contrast  to  this  is  a  narrow  region  at 
87  km  (<  1  km  widej  at  87  km  range  where 
near-zero  values  of  ZOR  extend  down  to  the 
ground.  We  believe  that  this  region  contains 
tumbling  hail  falling  to  the  ground 
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shown  in  Fig.  1-  In  Fig.  la,  contours  of  the  conventional  reflec¬ 
tivity  Z  are  drawn  evey  10  dBZ,  from  10  dBZ  to  70  dBZ  above 
1mm4  m'1.  In  Fig.  16,  contours  of  ZOR,  observed  simul¬ 
taneously.  are  drawn  at  1-dB  intervals.  The  dotted  Z  contour 
at  5  dBZ  in  a  and  6  delineates  the  overall  extent  of  the  echo 
On  this  day  the  0°C  isotherm  was  at  about  3  km  altitude,  and 
(as  is  typical  for  convective  clouds  )  values  of  ZDR  at  higher 
altitudes  in  Fig.  16  are  very  close  to  zero,  indicating  randomly 
tumbling  ice  particles.  Closer  to  the  ground  where  it  is  warmer 
than  0°C.  rain  is  to  be  expected,  and  as  the  ram  intensity 
increases  the  drops  become  larger  and  more  oblate,  leading  to 
larger  values  of  ZOR,  assuming  a  Marshall-  Palmer  raindrop  size 
distribution*,  calculations  show9  that  ZDR  should  exceed  1  dB 
when  the  Z  value  exceeds  35  dBZ  In  contrast  to  this  is  a  region 
at  a  range  of  about  87  km  where  low  values  of  ZDR  extend  down 
to  the  ground.  For  a  width  of  about  600  m  (that  is,  over  two 
neighbouring  samples),  a  ZDR  below  0.5  dB  accompanied  by  a 
Z  of  65  dBZ  extends  down  to  a  radar  elevation  of  zero  degrees. 
For  monodispersed  raindrops  this  value  of  ZOR  implies 
diameters  less  than  2  mm,  which  would  need  to  be  present  in 
the  unrealistic  concentration  of  500,000  to  give  the  observed 
Z  This  corresponds  to  rainfall  rates  of  thousands  of  millimetres 
per  hour,  and  any  other  distribution  would  require  even  higher 
concentrations.  The  only  plausible  explanation  is  that  the  pre¬ 
cipitation  is  hail  which  is  nearly  spherical  or  tumbles  as  it  falls. 
The  hail  must  be  large,  and  with  a  high  terminal  velocity,  so 
that  although  the  surface  is  wet  it  has  not  melted  sufficiently  to 
form  an  oblate-shaped  water  coating.  A  similar  but  less  extreme 
case  can  be  seen  in  Hall  et  ai 

To  monitor  hail  over  a  large  area,  the  radar  could  be  scanned 
in  azimuth  at  low  elevation  so  that  in  the  summer  only  precipita¬ 
tion  at  temperatures  much  above  0  9C  is  sampled.  If  regions 
where  Z  is  above  45  dBZ  are  considered,  then  any  region  where 
ZnR  ts  less  than  0.7  dB  can  be  identified  as  hail  reaching  the 
ground  Hatl-swaths  could  be  tracked  by  monitoring  the  move¬ 
ment.  m  successive  scans,  of  such  low-ZDR  cores  embedded  in 
regions  of  both  high  Z  and  ZDR. 

The  technique  outlined  above  would  not  result  in  false  alarms. 
However,  the  precise  shape  and  fall-mode  of  hail  is  not  well 
known  and  if  asphencal  hail  falls  with  a  degree  of  common 
alignment,  then  Z0R  might  have  a  local  minimum  but  not  fall 
below  the  0  7  dB  level  In  1986  the  Chilbolton  radar  will  be  able 
to  detea  the  orthogonal  return  in  addition  to  the  co-planar  one. 
This  information  may  provide  additional  corroborating  evidence 
for  hail.  The  linear  depolarization  ratio,  LOR,  for  vertically 
polarized  transmitted  radiation  is  defined  as 

TqR  =  10  logl0  i  Zv  h  Zv  i 


where  Zv  H  is  the  horizontally  polanzed  return  LOR  values  of 
less  than  -20dB  have  been  measured10  in  ram  at  a  3  2-cm 
wavelength.  The  orthogonal  return  is  very  low  in  amplitude 
because  it  is  believed11  that  raindrops  fail  with  their  axes  aligned 
in  the  vertical,  with  canting  angles  of  less  than  I*  Hail  is 
generally  asphencal  and  computations'^  show  that  randomly 
tumbling  wet  oblate  spheroids  with  an  axial  ratio  of  0.7  should 
have  an  LDR  greater  than  -20 dB  If  this  is  so.  the  CDR  hail 
signature  is  simple,  summertime  values  of  Ln„  above  -20  dB 
close  to  the  ground  must  be  due  to  hail.  The  extreme  case  of 
sphencal  hail  would  escape  LDR-detection,  but  would  give  a 
clear  zero  ZDR  signal  At  present  there  are  no  reports  of  measure¬ 
ments  of  Ldr  at  a  \0-cm  wavelength,  but  analysis  of  data  at 
3  cm  (supplied  by  P  H  Herzegh,  personal  communication > 
confirms  that,  although  propagation  problems  arfect  the  absolute 
values  of  LDR,  a  local  increase  of  about  6  dB  appears  to  coincide 
with  a  hail  shaft. 

The  presence  of  hail  seems  to  be  the  only  physically  realistic- 
explanation  of  the  localized  region  of  low  ZDR  embedded  in  a 
more  extensive  region  where  both  Z  and  Zt)n  are  high  Although 
it  has  not  been  possible  to  obtain  independent  ground-based 
observations  in  the  500-m-wide  region  shown  in  Fig  1.  a  long 
senes  of  comparisons13  of  the  radar  data  with  simultaneous 
raindrop  size  distributions  measured  at  the  ground  reveals  that 
such  high  Z-low  ZDR  combinations  are  not  found  in  rainfall 
More  extensive  observations  with  the  Chilbolton  radar  scanning 
in  azimuth  at  low  elevation  should  enable  the  tracking  of  low- 
ZDR  cores  and  any  accompanying  anomalous- regions,  and 
thus  increase  the  probability  of  obtaining  independent  ground 
based  confirmation  of  the  hail  identification 
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l.  introduction 

As  raindrops  grow  in  size  they  become  more 
oblate  and  consequently  the  radar  reflectivity  factor 
for  horizontally  polarised  radiation  (2H)  becomes 
increasingly  larger  than  that  for  the  vertical.  The 
differential  reflectivity,  ZDR,  expressed  in  decibels  is 
defined  as 

2 OR  =  10  log  (ZH/ZV) 

ZDR  is  positive  for  ram  and  its  magnitude  provides  an 
estimate  of  th*  mean  size  of  the  raindrops.  Because 
ZDR  is  independent  of  concentration,  observations  of 
Z  and  ZDR  can  be  fitted  to  any  two  parameter  rain¬ 
drop  size  distribution  and  used  to  derive  more  accurate 
rainfall  rates  than  those  obtained  from  2  alone  (Seliga 
and  8ringi,  1976,  Goddard  and  Cherry,  1904). 

The  Chilbolton  radar,  etuated  in  Hampshire 
England,  operates  at  9.75  cm  with  a  quarter  degree 
beamwidth  and  is  able  to  measure  Z  to  within  0.7d8 
and  in  widespread  precipitation  ZDR  of  0.1d8  (Cherry 
and  Goddard,  1982).  If  we  consider  raindrops  as  oblate 
spheroids  and  apply  Cans'  scattering  theory,  then  a 
change  of  ZOR  of  Q.WB  is  equivalent  to  a  change  in 
axial  ratio  of  less  than  0.01.  The  best  currently 
available  laboratory  measurements  of  equilibrium 
raindrop  shapes  (Pruppacher  and  Pitter,  1971)  have 
error*  in  axial  retie  considerably  greater  than  0.01.  In 
a  seriea  of  comparison*  of  radar  observations  for 
values  of  ZDR  up  to  2d8,  with  both  ground  based 
disdrometers  (Goddard  et  at,  1983)  and  airborne 
measurements  (Cherry  et  al,  1984)  found  significantly 
improved  agreement  if  they  used  slightly  more 
spherical  shapes  than  those  of  Pruppacher  and  Pitter. 
The  changes  proposed  by  Goddard  and  Cherry  are 
actually  within  the  experimental  error  of  the 
Pruppacher  and  Pitter  measurements. 

A  ZDR  value  of  2d8  is  equivalent  to  a  mean  drop 
size  of  about  1.7mm,  and  average  rainfal'  rates  of 
about  40mm/hr.  In  this  paper  we  analyse  the  values  of 
ZDR  for  heavier  rainfall  rates  and,  so  that  consistency 
may  be  preserved  with  the  radar  results,  suqgest  an 
empirical  modification  to  the  Pruppacher  and  Pitter 
shapes  for  the  larger  drop*.  In  the  light  of  these 
modifications,  we  then  interpret  some  observations  of 
the  evolution  of  ZOR  within  vigorous  convective 
clouds  in  terms  of  the  size  and  concentrations  of  rain¬ 
drops.  We  have  payed  particular  attention  to4the 
narrow  columns  of  positive  2DR  which  occasionally 
extend  above  the  zero  degree  isotherm  during  vigorous 
convection. 


2.  THE  SHAPE  OF  LARGE  RAINDROPS. 

Ooviak  and  Zmie  (1984)  have  calculated  that 
ZDR  values  above  3.5dB  imply  targets  other  than  rain¬ 
drops  and  suggest  that  they  result  from  ice  particles 
which  have  very  oblate  shape*  when  they  melt.  Thu 
assertion  is  based  on  Green's  (1975)  formula  for  the 
axial  ratio  of  raindrops,  which  is  very  close  to  the 
measurements  of  Pruppacher  and  Pitter.  Using  these 
shapes  a  6mm  drop  would  be  associated  with  a  ZDR  of 
3.8d8  and  an  0mm  drop  4.7d8;  summing  the  weighted 
contributions  of  the  various  drops  in  a  Marshall-Palmer 
raindrop  size  distribution  leads  to  the  conclusion  that 
the  ZDR  of  rain  cannot  exceed  3.8d8. 

Values  of  ZDR  above  3.5dB  and  as  high  as  6dB 
have  been  reported  in  convective  showers  (Illingworth 
et  al,  1984)  and  in  view  of  the  above  argument  it 
seems  logical  to  attribute  such  observation  to  meltinq 
ice  particles.  On  closer  examination  other  possibilities 
arise.  Pint,  recent  laboratory  measurement*  of  melt¬ 
ing  particles  suggest  that  unless  the  ice  particles  are 
initially  extremely  oblate  then  the  shapes  during 
melting  are  no  more  distorted  than  those  of  the  same 
sized  raindrop,  and  second,  there  are  considerable 
experimental  errors  in  the  Pruppacher  and  Pitter 
shapes  for  the  larger  drops. 

The  laboratory  studies  of  Pasmussen  et  al  (1984) 
show  that  as  spherical  hail  particles  melt  they  become 
progressively  more  oblate  until,  on  complete  meltinq, 
their  axial  ratio  la  equal  to  that  of  the  equivalent 
water  drop.  Hail  pellets  larger  than  9mm  always  shed 
the  excess  water.  A  ZOR  of  6d8  could  result  from 
horizontally  alityied  ice  particles  with  a  wet  surface 
and  an  axial  ratio  of  0.52.  This  seems  unlikely  because 
such  large  solid  particles  should  tumble,  and  hail  is 
rarely  observed  to  have  such  an  oblate  shape.  High 
values  of  ZDR  attributed  to  melting  snowflakes  (Hall 
et  al,  1984)  would  only  occur  in  stratiform  clouds  and 
should  be  restricted  to  the  bni^t  band. 

There  is  considerable  uncertainty  as  to  the 
precise  shape  of  raindrop*  above  4mm  diameter. 
Laboratory  experiments  are  difficult  because  the  drops 
oscillate,  and  the  degree  of  oscillation  is  probably 
dependent  upon  the  amount  of  turbulence  in  the  wind 
tunnel.  Pruppicher  and  Pitter  quote  equilibrium  axial 
rptios  for  drops  of  6  and  8mm  diameter  which  would 
lead  to  ZDR  values  of  3.8  and  4,7d8  respectively,  but, 
because  o'  the  experimental  errors  quoted,  these  axial 


ratio*  could  be  as  low  as  0-62  and  0.51,  implying  a  ZDR 
of  4.7  and  6.5d8.  It  is  difficult  to  check  these  pre¬ 
diction*  m  the  field.  Direct  aircraft  measurements  by 
Cooper  et  aJ  (1983)  showed  that  for  a  given  equivalent 
diameter  there  was  a  spread  of  axial  ratios  implying 
appreciable  oscillation,  but  that  the  mewi  shape  for 
5mm  rt ops  was  significantly  more  distorted  than 
°ruppacher  and  PitteCs  measurements.  Only  «  yerv 
few  drops  above  6mm  were  sampled  but  they  indicated 
an  axial  ratio  of  onlv  0.5  for  a  6.3mm  drop. 

If  the  concentration  of  raindrops,  N,  of  diameter 
O,  is  assumed  to  follow  the  Marshall-Palmer  raindrop 
distribution 

N(D)  *  No  exp  (-3.67  Oo/D) 

where  No  is  8000/m  J/mm,  then  if  the  drop  shapes  are 
known  *e  can  calculate  ZDR  as  a  function  0f  Z. 
figure  1  compares  the  predicted  variation  using  the 
Pruppacher  and  Pitter  shapes  and  an  extrapolation  of 
the  more  oblate  measurements  of  Cooper  et  al.  The 
actual  slues  of  axial  ratio  used  in  these  computations 
are  displayed  in  Table  1.  It  should  be  stressed  that 
Cooper  et  al  did  not  measure  any  drops  above  7mm, 
and  that  the  axial  ratios  for  drops  above  this  size  are 
obtained  by  linear  extrapolation.  The  dotted  curves  in 
Figure  !  were  obtained  by  truncating  the  exponential 
distribution  at  3mm;  observation  shows  wi  icant 

number  of  drops  larger  than  this  size.  Such  a  trunca¬ 
tion  can  be  introduced  more  naturally  bv  usmq  a 
gamma  distribution  fUlbrich,  1983),  and  tr«  dashed 
branches  of  the  curves  are  the  predicted  variations  for 
an  m  of  2  in  the  gamma  distribution.  The  solid  branch 
of  the  curves  is  obtained  bv  truncating  a  pure  exponen¬ 
tial  (m  j  Q)  at  7mm.  These  variations  introduced  bv 
different  choices  of  drop  size  distribution  are  small 
compared  with  the  effect  of  the  two  different  models 
for  the  ,*op  shapes.  If  the  more  oblate  shapes  are 
adopted  then  it  appears  that  values  of  ZDR  of  up  to 
6dB  can  arise  in  heavy  ram. 


3.  RADAR  (XRJVED  SHAPES  OF  IARC£  DROPS. 

We  have  indicated  above  that  independent,  veri¬ 
fication  of  the  predicated  Z/ZDR  relationship  by 
direct  measurement  of  the  raindrops  within  the  radar 
beam  is  not  practically  possible.  Iratead  we  have 
adopted  a  statistical  approach;  althomji  on  individual 
occasions  the  crop  size  distributions  we  not  M*z-jhall- 
Palmer,  if  sufficient  data  points  are  taken  then  the 
average  variation  of  Z  as  s  function  of  ZDR  sTviuld 
provide  an  indication  of  the  drop  shapes.  Over  ninety 
thousand  data  points  from  storms  observed  on  six 
different  days  m  1985  are  included  m  the  plot  dis¬ 
played  m  figure  2,  where  the  average  values  of  ZDR 
have  been  calculated  for  every  2d8  increase  in  Z.  for 
comparison  purposes  the  theoretical  curves  from 
f  igure  1  for  the  two  different  sets  of  drops  shapes 
predicted  for  exponentials  truncated  at  8mm  are  also 
included.  On  these  days  the  freezing  level  was  at  3km 
and  all  data  points  for  altitudes  between  1.25  and  1km 
were  analysed.  At  such  low  altitudes  most  of  the 
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fable  1:  The  axial  ratio  (b/ a)  and  differential 

reflectivity  (ZDR)  for  raindrop*  having  equivalent 
diameter  X))  of  4  to  8mm  using  the  drop  shapes  of 
Pruppacher  and  Pitter  (1971'  and  values  extrapolated 
from  Cooper  et  al  (1983). 
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4.  EVOLUTION  OF  RAINDROPS. 


precipitation  should  be  rain,  but  on  a  few  scans  the 
Z/ZDR  plots  showed  a  fall  in  ZDR  for  values  of  Z 
above  50dBZ.  Figure  3  shows  ar  example  of  this.  We 
believe  {Illingworth  e<  si,  1986/  that  this  behaviour  is 
due  to  large  tumblirxj  hail  particles  surviving  the  2Vm 
fall  below  the  freezing  level.  Figure  4  i*  a  contour 
plot  of  a  PPI,  clearly  showing  the  core  of  ZDR  below 
Q_5dB  occurmg  at  tb*  centre  of  the  maximum  Z  echo, 
which  we  interpret  as  a  had  shaft.  The  occasional 
section  showing  such  characteristics  was  excluded 
from  the  analysis.  The  standard  deviation  0f  the 
points  in  Figure  2  for  Z OR  above  2dB  is  about  0.7  d8  in 
ZDR,  showing  that  on  the  average  each  mdiviAjal  size 
distribution  does  not  have  a  value  No  s  8000. 

T>w  observed  variation  of  ZDR  with  Z  in  natural 
ram  shown  in  Figure  2  deviates  significantly  from  the 
predictions  of  the  Pruppacher  and  Pitter  shapes,  and  ts 
very  dose  to  the  values  computed  using  the  much 
more  oblate  shapes  extrapolated  from  Cooper  et  al. 
An  alternative  explanation  would  be  to  consider  that 
the  Pruppacher  and  Pitter  shapes  are  valyi,  but,  that 
when  Z  it  high.  No  falls  below  8000/m  /mm.  For 
example  a  ZDR  of  3.6d8  should  on  average  be 
accompanied  by  a  Z  of  70d8Z  if  No  is  8000/m  / mm.  In 
fact  the  average  value  ;S  only  58d8Z  implying  that,  for 
this  high  value  of  Z,  No  is  less  than  800,  and  that  the 
rainfall  rate  is  correspondingly  lower  than  that  which 
would  occur  for  No  *  3000.  This  seems  unlikely. 
Calculated  rainfall  rates  using  No  =  8000  agree  with 
the  empirical  Z/R  relationships  to  within  !d8.  If  we 
use  the  more  oblate  ramdroo  shapes  in  the  model  then 
the  empirical  Z/R  relationahip  is  not  disturbed,  the 
observed  variation  of  Z  with  ZDR  can  be  explained, 
and  furthermore  the  high  values  of  ZDR  of  6dB  can  be 
attributed  to  raindrops  without  invoking  extraordinary 
oblate  aligned  melting  ice  particles. 


one  PPr.  ^uch  data  is  attributed  to  hail  and  is 
not  included  in  Fi®  2.  dotted  lines  as  in  Fin  2. 


In  convective  clouds  the  value  of  ZDR  above  the 
0°C  isotherm  is  observed  to  be  very  clcae  to  zero  with 
values  generally  below  0. ldB,  and  this  is  interpreted  as 
tumbling  grsupel.  On  occasion  localised  ereas  of 
positive  ZDR  extend  above  the  freezing  level 
"Illingworth  et  al,  1904)  and  these  are  generally  in  the 
form  of  narrow  column#  as  shown  m  Figure  5. 

During  tr»  summer  of  1903  twenty  columns  of 
♦  ve  ZDR  having  values  of  ZDR  above  2d8  were 
observed.  In  ten  of  these  cases  it  was  possible  to  trace 
the  evolution  of  the  column  from  three  or  more 
sequential  RHTs  taken  at  the  same  azimuth,  each  RHI 
et  the  same  axjmuth  being  typically  separated  by 
about  two  minutes.  These  esses  occur ed  in  both 
isolated  initial  echoes  and  embedded  within  mature 
cumulonimubs  clouds.  Considering  the  individual  RHJa, 
we  conclude  that  the  positive  columna  of  ZDR  have 
the  following  characteristics: 

a)  The  maximum  height  of  the  2dB  ZDR  contour  was 
1.7  km  above  the  zero  degree  isotherm  where  the 
temperature  was  about  -10^C.  Just  above  this 
heicftt  m  the  reqion  su*rondinq  the  column 
ZDR  values  are  less  that  0.  ldB. 

b)  Well-developed  columna  rising  more  than  700m 
above  the  freezing  level  have  a  circular  cross- 
section  with  a  diameter  bet  ween  1  and  2Vm. 

c)  F or  such  columna  the  maximum  value  of  ZDR 
exceeds  4.5dB  and  can  reach  6dB;  a  further  500m 
higher,  values  of  ZDR  are  zero. 


observed  at  1517  on  8/Jul/83.  The  hatched  area 
where  Z  is  below  h.5dB  and  2  is  above  85dB7  is 
attributed  to  hail.  Contours  in  ZD*  are  every  ldB 
startinc  at  3.5 d8  and  in  2  are  ev*rr  IhdB. 
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5. 


CONCLUSIONS 


From  the  sequential  RHIs  we  conclude: 

d)  A  well  defined  positive  column  doe*  not  persist 
for  more  than  ten  minutes. 

e)  While  the  positive  column  ts  present  the 
magnitude  of  2  increases  try  up  to  20dBZ,  finally 
exceeding  58d8Z.  In  one  case  however  the  value 
of  2  only  rose  from  >8  to  66dBZ. 

0  During  the  development  of  the  positive  column 
the  *2dB  2 DR  contour  rises  at  about  7m/s;  this 
rise  being  coincident  with  a  similar  rate  of  ascent 
for  either  the  20dBZ  contour  in  2  for  the  isolated 
developing  cells,  or  the  40d8Z  contour  for  the 
columns  embedded  in  existing  ceils.  These  figures 
imply  that  the  air  itself  is  rising  at  about  15m/s. 

g)  Interpreting  the  Z  and  ZDR  data  in  the  column  as 
the  size  and  concentration  of  supercooled 
raindrops,  the  derived  values  of  No  always  tended 
to  increase  with  time  towards  Marshall- 
Palmer  value  of  8000.  For  isolated  convection 
this  increase  in  6  minutes  was  about  a  factor  of 
50,  and  for  the  embedded  convection  a  factor  of 
two  in  two  minutes.  In  the  one  case  (see  (e))  No 
remained  near  200. 

h)  After  the  disappearance  of  the  positive  ZDR 
column,  the  position  of  the  maxima  in  both  Z  and 
ZDR  fell  in  altitude  at  a  velocity  consistent  with 
the  terminal  velocity  of  the  large  precipitation 
particles. 


Fig  5:  An  WI  on  6/Jtil/*3  at  1651  showing  a 
colurtn  of  ♦v*  2DP  reaching  a  height  of  i.TYre. 

P  C  isochens  3.2lua.  Contours  in  Z  every  lOdBZ 
starting  at  5dB2,  and  in  2DP  every  1.9dB  starting 
at  0.5dB.  Hatched  area  Z  above  55dBZ,  ZDR  above 
4 . 5dB. 


The  observed  values  of  Z  and  ZDR  m  fwavy 
precipitation  are  consistent  with  the  more  oblate 
shapes  reported  by  Cooper  et  ei.  We  suggest  that 
positive  value*  of  ZDR  of  up  to  6dB  are  due  to  large 
raindrops  which  may  be  supercooled,  and  that  hail  or 
qraupel  particles  which  are  wet  will  have  smaller 
values  of  ZDR  than  the  equivalently  sized  raindrop. 
Transient  positive  columns  of  ZDR  extending  up  to 
1,7km  above  the  freezing  level  occur  dunnq  vigorous 
convection;  initially,  raindrop  concentrations  are  low 
but  they  evolve  towards  the  Marshall-Palmer  distribut¬ 
ion. 
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APPENDIX  C 


INTRODUCTION 


Our  present  knowledge  of  how  precipitation  develops  within  clouds  has  been  obtained 
using  both  radar  and  instrumented  aircraft.  Measurements  of  the  radar  reflectivitv, 
Z,  provide  an  indication  of  the  rainfall  rate;  air  motions  can  be  computed  from  multiple 
doppler  radar,  and  hydrometeor  types  and  sizes  may  be  obtained  from  probes  mounted 
on  aircraft.  However,  information  on  hydrometeors  is  generally  available  only  for  a 
very  limited  number  of  one  dimensional  penetrations  through  the  less  vigorous  regions 
of  a  cloud,  and  a  complete  history  of  the  hydrometeor  evolution  throughout  the  volume 
of  the  cloud  can  only  be  built  up  by  extrapolation  from  these  limited  samples  (e.g.  Hallett 
et  al,  1978;  Knight  and  Squires,  1982;  Heymsfield  and  Hjelmfelt,  1984,  and  for  hail  growth. 
Browning  et  al,  1976,  Browning  1978). 

Several  radar  scans  through  a  cloud  do  provide  a  reasonably  complete  three  dimensional 
picture  of  the  extent  and  intensity  of  the  radar  reflectivity,  2,  but  2  itself  gives  no 
information  on  hydrometeor  types,  sizes  or  concentrations.  Several  empirical  relation¬ 
ships  are  available  which  relate  2  to  the  rainfall  rate  R,  but  because  2  is  proportional 
to  the  product  ND^  summed  over  particles  of  all  sizes  (where  N  is  the  concentration 
of  hydrometeors  of  diameter  D),  the  use  of  such  2-R  relationships  is  error  prone  and 
is  equivalent  to  assuming  that  all  hydrometeor  size  distributions  are  identical.  The 
same  2  could  result  from  many  different  size  distributions  having  widely  varying  equiva¬ 
lent  rainfall  rates.  An  additional  ambiguity  occurs  if  mixed  phase  clouds  are  present, 
because  liquid  water  has  a  reflectivity  6.7dB  greater  than  ice;  for  melting  particles 
interpretation  of  2  is  even  more  error  prone. 


In  this  paper  we  present  measurements  not  only  of  2  but  of  the  differential  radar  reflecti¬ 
vity,  of  evolving  convective  clouds.  This  additional  parameter  can  reduce  consider¬ 

ably  the  ambiguities  present  in  the  interpretation  of  Z  alone. 
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measurement.  Each  range-gate  data  sample  comprises  a  spatial  average  (linear 
power)  over  300m  (four  pulse  volumes).  At  60km  the  sample  volume  is  typically 
a  cube  of  side  300m  in  length.  Time  averaging  (linear  power)  for  each  data  sample 
is  over  210ms  every  240ms  (i.e.  64  transmitted  pulses  on  each  polarisation).  The 
worst  case  systematic  errors  in  measuring  Z  and  Zpjp  were  estimated  to  be  O.TdB 
and  O.ldB,  respectively,  while  the  standard  deviations  of  the  random  errors  when 
measuring  rain  echoes  were  estimated  to  have  the  same  values.  A  detailed  study 
of  the  accuracy  of  the  measurement  has  been  given  by  Cherry  and  Goddard  (1982). 

Herzegh  and  Carbone  (1984)  have  shown  that  when  strong  reflectivity  gradients 
in  Z  are  present,  a  mismatch  of  the  polarisation  characteristics  of  the  antenna 
sidelobes  may  give  rise  to  totally  false  Zqr  values.  In  this  case  it  is  possible  for 
the  main  beam  to  be  sensing  a  low  reflectivity  region,  but  for  a  significant  part 
of  the  reflected  signal  to  be  due  to  the  sidelobe  detecting  the  high  Z  region.  They 
give  examples  from  the  1°  beamwidth  CP2  radar  where  mismatched  sidelobes  result 
in  spurious  regions  of  high  positive  and  negative  Zqr  on  either  side  of  a  reflectivity 
maximum.  Very  much  higher  reflectivity  gradients  would  be  required  to  produce 
such  effects  for  the  quarter  degree  beamwidth  of  the  Chilbolton  radar,  but  we 
shall  return  to  this  potential  problem  when  discussing  the  observations  in  Section  3. 

(b)  Zqr  of  Rain 

If  hydrometeors  are  assumed  to  be  oblate  spheroids,  having  their  figure  axes  vertical, 
then  application  of  -'■ns  scattering  theory  gives  the  curves  in  Figure  1  (from  Hall 
et  al,  1984)  showing  .ie  value  of  Zqr  as  a  function  of  axial  ratio  for  liquid  water 
and  air-ice  mixtures.  The  Figure  clearly  demonstrates  that  for  a  given  shape  this 
ratio  is  much  larger  for  liquid  water  than  for  ice.  For  rain  it  is  probably  valid  to 
assume  that  the  drops  fall  with  their  figure  axes  aligned  in  the  vertical,  but  any 
canting  or  tumbling  will  lead  to  lower  values  of  Zqr:  in  the  extreme  case  of  random 
tumbling,  which  may  occur  for  ice  particles,  Zpp  will  be  zero. 
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The  most  accurate  laboratory  measurement  of  equilibrium  raindrop  shapes  is  provided 
by  Pruppacher  and  Pitter  (1971).  In  a  series  of  comparisons  of  radar  observations 
for  values  of  Zqr  up  to  2dB,  with  both  groundbased  distrometers  (Goddard  et  al, 
1983)  and  airborne  measurements  (Cherry  et  al,  1984)  Goddard  and  Cherry  (1984a) 
found  significantly  improved  agreement  if  they  used  slightly  more  spherical  shapes 
than  those  of  Pruppacher  and  Pitter.  This  modified  relationship  between  drop  diame¬ 
ter  and  Zpp  is  displayed  on  the  abscissa  of  Figure  2.  Oscillations  of  raindrops 
are  believed  to  be  asymmetric,  and  if  these  are  excited  then  a  time  averaged  shape 
should  be  more  spherical  than  the  equilibrium  one,  and  this  could  reduce  the  observed 
value  of  Zqr  (Beard  et  al,  1983).  It  could  be  argued  that  because  oscillations  are 
collision  induced,  any  correction  to  Zqr  would  have  to  be  concentration  dependent. 
Direct  photography  of  raindrops  at  the  ground  (Jones,  1959)  indicates  that  a  degree 
of  canting  is  present,  but  the  wind  shear  at  the  ground  may  not  be  representative 
of  conditions  within  clouds;  computations  (Beard  and  Johnson,  1983)  suggest  that 
shear  produced  by  turbulence  would  lead  to  small  canting  angles  distributed  around 
a  mean  zero  value.  However,  the  changes  proposed  by  Goddard  and  Cherry  are 
actually  within  the  experimental  error  of  the  Pruppacher  and  Pitter  measurements, 
and  so  there  is  no  evidence  that  canting  or  collision  induced  oscillations  affect 
the  Zqp  measurement  of  rain.  The  calibration  of  Zpjp  is  so  sensitive  to  drop  shape 
that  at  present  it  is  able  to  detect  changes  in  shape  that  cannot  be  resolved  by 
direct  photography. 

If  both  Z  and  Zqr  are  measured  for  rain  then  the  data  can  be  fitted  to  any  two 
parameter  raindrop  distribution.  Figure  2  considers  one  extreme  example,  the 
monodispersed  spectrum,  and  for  a  given  Zqr  the  curves  show  the  expected  values 
of  Z  for  a  rainfall  rate  of  1mm  hr"l  and  for  drop  concentrations  of  1  and  100  per 
cubic  meter.  Because  Zgp  is  independent  of  concentration,  the  values  of  N  and 
R  scale  linearly  with  Z.  There  is  uncertainty  over  the  average  shape  of  drops  having 
an  equivalent  diameter  above  4mm.  Some  recent  aircraft  measurements  (Cooper 
et  al,  1983)  suggest  that  they  are  more  oblate  than  the  equilibrium  shapes  of 
Pruppacher  and  Pitter,  with  a  6mm  drop  having  a  Zpjp  of  5.8dB  instead  of  3.8dB, 


and  an  indication  that  the  values  of  7mm  drops  should  be  7.7dB  rather  than  4.18dB. 
The  two  dotted  branch  to  the  curves  in  Figure  2  for  Zqr  above  2.5dB  are  computed 
using  these  more  oblate  shapes. 

The  Marshall  Palmer  raindrop  size  distribution  is  given  by 

N(D)  =  N0  (exp(-3.67  D/D0))  (2) 

where  N0  is  8000m'^  mm'l  and  constant,  and  D0  depends  upon  the  rain  rate  and 
is  equivalent  to  a  median  drop  size  such  that  an  equal  volume  is  contained  in  drops 
smaller  and  larger  than  DQ.  If  both  Z  and  Zq^  are  available  then  a  fit  to  a  general¬ 
ised  exponential  distribution  with  both  N0  and  D0  as  variables  is  possible  (Seliga 
and  Bringi,  1976).  In  Figure  3  the  abscissa  shows  the  relationship  between  Zqr 
and  D0  obtained  by  weighting  the  Z-D  dependence  in  Figure  2  and  truncating  the 
exponential  spectra  at  8mm  diameter.  As  above,  because  they  scale  linearly  with 
Z,  the  values  of  N0  and  R  may  be  obtained  by  linear  extrapolation  of  the  plotted 
curves.  According  to  this  Figure  values  of  Zqr  above  5dB  cannot  be  fitted  to 
realistic  exponential  raindrop  size  distributions  if  the  drop  shapes  of  Pruppacher 
and  Pit  ter  are  used. 

Some  caution  should  be  exercised  in  interpreting  these  curves  for  Zpp  values  above 
3dB  because  of  the  uncertainty  in  the  precise  shapes  of  the  larger  drops  and  the 
sensitivity  of  the  curves  to  the  truncation  of  the  exponential.  If  the  exponential 
is  not  truncated  at  8mm  then,  for  D0  above  4mm,  because  the  very  few  drops  above 
8mm  are  so  distorted,  their  contribution  determines  Zp>£.  Consequently  in  very 
intense  rain  the  shape  of  the  drops  at  the  truncation  limit  of  8mm-fixes  the 
maximum  possible  value  of  Zqr  to  about  3.5dB  for  Pruppacher-Pitter  shapes. 
If  the  more  oblate  shapes  of  Cooper  et  al  are  used  then  an  8mm  truncation  limit 
would  lead  to  a  maximum  Zqr  of  around  6dB. 
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The  sensitivity  of  predicted  values  of  Zj^p  in  very  intense  rain  to  the  rather  arbitrary 
choice  of  the  truncation  limit  of  the  exponential  is  rather  artificial;  the  observed 
concentrations  of  these  very  large  drops  are  in  fact  smaller.  Ulbrich  and  Atlas 
(1984)  have  suggested  introducing  this  effect  more  naturally  by  the  use  of  gamma 
distributions  of  the  form: 

N(D)  =  N0Dm  exp(-<3.67  +  m)  D/D0})  (3) 

The  third  variable,  m,  is  zero  for  a  pure  exponential  and  infinity  for  monodispersed 
spectra  and  is  a  measure  of  spectral  width.  Ulbrich  and  Altas  (1984)  suggested 
using  m  =  2,  but  Goddard  and  Cherry  (1984b)  obtained  even  better  agreement  with 
measured  rainrates  for  a  value  of  m  =  5.  The  magnitudes  of  N,  N0  and  R  calculated 
from  Figures  2  and  3  (m  =  infinity  and  zero)  represent  limits,  with  the  true  value 
of  the  total  drop  concentration  probably  lying  between  these  two  extremes.  We 
note  that  for  an  observed  Z  and  Zqr  combination  the  choice  of  the  monodispersed 
spectrum  will  result  in  predicted  rainfall  rates  typically  3dB  greater  than  for  the 
exponential  spectrum. 

(c)  Measurements  of  Hydrometeors  above  the  0°  isotherm 

The  shape  of  an  ice  particle  and  its  mode  of  fall  is  not  a  unique  function  of  its 
size  and  so  the  analysis  of  the  signal  is  complicated.  Hall  et  al  (1984)  consider 

the  problem  in  detail,  and  Figure  1  taken  from  their  paper,  demonstrates  some 
of  the  possible  ambiguities.  Even  if  the  particles  were  ail  aligned  with  the  major 
axis  horizontal  we  see  that  Zqr  would  depend  not  only  on  the  particle  shape  but 
also  on  the  density.  If  they  are  not  perfectly  aligned  as  they  fall,  then  the  value 
of  Zqp  will  be  lower,  falling  to  zero  for  random  tumbling.  Further  complications 
and  ambiguities  would  arise  if  the  actual  ice  particle  shapes  were  used  rather  than 
idealised  oblate  spheroids. 


From  experience  with  the  Chilbolton  data  Zqr  is  usually  very  ■  lus>-  \  '  z.-rn  ,r, 
regions  of  convective  clouds  where  the  temperature  is  below  zero,  mill;  .itmt  m*  ,rlv 
spherical  or  randomly  orientated  tumbling  ice  particles.  In  the  next  tior.  we 
shall  be  considering  transitorv  regions  of  high  positive  Zqr  above  4dB  within  con¬ 
vective  clouds  which  extend  both  above  and  beiow  the  zero  degree  isotherm 
occur  where  Z  is  moderate  or  high.  We  shall  now  consider  the  tvpes  of  hvdrometeor 
which  could  give  rise  to  such  signals. 

(i)  Dry  ice  particles 

From  Figure  1  snow  cannot  give  a  Zqr  above  ldB,  but  a  Zqr  of  -4dB  can 
result  from  dry  aligned  oblate  spheroids  of  ice  having  axial  ratios  of  0.15 
for  an  ice-air  mixture  of  demity  0.5g  cm'^  or  spheroids  of  solid  ice  with 
a  ratio  of  0.4.  For  example.  Hall  et  al  (19841  found  regions  of  high  positive 
Zqr  in  the  upper  regions  of  stratiform  clouds  where  Z  is  low  and  suggested 
that  these  may  be  due  to  small  ice  crystals  in  the  form  of  plates  which  fall 
with  their  larger  dimension  horizontally  aligned.  Direct  aircraft  observations 
(Bader  et  al,  19861  flying  through  such  clouds  in  1984  have  confirmed  this 
suggestion.  However,  such  small  particles  can  give  only  low  values  of  Z  and 
on  melting  would  acquire  a  water  coating  and  give  much  larger  values  of 
Zqr  just  below  the  0°C  isotherm  before  subsequently  collapsing  to  smaller 
raindrops.  Hall  et  al  11984)  do  find  regions  of  positive  Zqr  in  association 
with  the  bright  band  m  Z  which  are  obviously  caused  by  melting  particles, 
but"  for  the  observations  we  shall  present  there  was  no  bright  band  and  the 
positive  Zqr  extended  2km  above  the  zero  degree  isotherm. 

The  various  habits  of  ice  crystal  growth  can  result  in  extremely  oblate  shapes 
but  because  of  their  size  would  be  associated  with  low  values  of  Z.  Larger 
ice  particles  are  formed  by  aggregation  or  riming  and  there  is  no  evidence 
from  aircraft  measurement  that  they  would  have  both  such  a  degree  of  oblate¬ 
ness  and  also  be  aligned  when  they  fall.  Although  the  fall  mode  of  larger 
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frozen  hydrometeors  and  hailstones  is  not  well  known  (Thwaites  et  al,  1977, 
Matson  and  Huggins,  1980),  most  evidence  suggests  that  discs  tumble  when 
the  Reynolds  number  exceeds  1000  (Pruppacher  and  Klett,  1980). 

(ii)  Ice  in  wet  growth 

Hail  or  graupel  particles  in  wet  growth  will  have  a  Zqr  of  4dB  if  the  axial 
ratio  is  0.65.  This  could  happen  in  two  ways;  either  sufficient  water  builds 
up  on  a  spherical  particle,  or  a  surface  layer  of  water  on  a  quite  asymmetric 
hail  pellet  stabilises  it  and  stops  it  tumbling.  The  particles  would  probably 
shed  any  excess  water,  but  if  for  the  sake  of  argument,  we  extend  the  analysis 
of  Bailey  and  Macklin  (1968)  to  smaller  particles  with  a  fall  velocity  of  V 
=  1.43  (D)0‘®  where  D  is  in  mm  and  V  in  m  s‘  1  (Pruppacher  and  Klett,  19801, 
we  find  that  if  the  surface  is  to  be  wet,  then  for  diameters  of  2,  4  and  10mm, 
a  cloud  liquid  water  content  of  13,  7.3  and  3.2g  m"^  respectively  is  required. 
Even  higher  values  would  be  needed  for  a  reasonable  thickness  of  water  to 
build  up,  and  for  rough  particles  with  higher  ventilation  coefficients  these 
values  should  be  increased  threefold.  From  sonde  ascents  made  on  the  davs 
in  question  the  maximum  adiabatic  liquid  water  content  at  0°C  and  -10‘ C 
cannot  exceed  5g  and  7g  m  3f  We  conclude  that  small  graupel  particles  in 
wet  growth  require  unrealistic  liquid  water  contents  and  would  also  melt 
to  give  much  smaller  Zqr  below  the  freezing  level.  Computation  (Mason, 
1956)  shows  that  for  solid  ice  spheres  to  fall  2km  below  the  freezing  level 
without  melting  completely,  the  original  diameter  must  be  at  least  4mm; 
for  graupel  of  density  0.3g  cm'^  the  equivalent  size  is  8mm.  However, 
interpretation  of  these  signals  in  terms  of  hail  of  diameter  greater  than  5mm 
which  is  in  wet  growth  cannot  be  ruled  out. 


(iii I  Supercooled  raindrops 


A  simpler  explanation  of  such  signals  would  be  to  assume  they  are  due  to 
large  supercooled  raindrops,  a  Zpp  *■  4dB  being  equivalent  to  a  diameter 
of  about  6mm;  and  when  analysing  these  results  we  shall  adopt  this  approach. 
We  shall  return  to  consider  the  wet  growth  aspects  in  the  discussion  Section. 

3  RESULTS  OBTAINED  IN  1983 

The  aim  of  the  1983  studies  was  to  trace  the  development  of  isolated  convective  clouds. 
Tvpical  echoes  from  these  clouds  have  horizontal  diameters  of  about  5km  and  about 
5  vertical  sections  (RHIs)  at  different  azimuths  were  sufficient  to  identifv  the  maximum 
echo  and  its  spatial  extent.  Each  RHI  takes  about  30  seconds,  so  an  observation  at 
any  particular  azimuth  was  obtained  every  Z  or  3  minutes.  Occasional  horizontal  sections 
(PPIs)  served  to  confirm  the  isolated  character  of  the  storms.  Hall  et  al  (1980!  reported 
that  pillars  of  -veZ^p  extending  higher  than  the  0:C  isotherm  were  occasionally  observ¬ 
ed,  and  suggested  that  thev  were  due  to  supercooled  raindrops.  One  particular  goal 
was  to  find  out  if  these  columns  of  large  -veZj^p  were  a  common  occurrence  and  to 
study  their  persistence.  Because  there  is  no  area  ol  preferred  convection  over  Southern 
England  the  location  of  first  echoes  is  difficult,  but  in  these  examples  the  problem  was 
alleviated,  because  winds  were  very  light  leading  to  minimal  advection  of  cells,  and 
on  6  anti  7  July  a  slight  shear  resulted  m  the  'development  ot  sequences  of  daughtei 
cells. 


3  l  a  u  y  .  6  i  r  o  6  5  5 

Or:  this  day  the  1000  hrs  stride  ascent  f  r  -  Craw  .  t  y , 
east  rf  Chi  Ido  1  ton,  showed  mat  a  further  r  ;se  surface 

temperature  of  0  or  i  decrees  would  o  e  sufficient  t  c  lift  t  r.  e 
moist  low  level  air  to  the  condensation  level  of  l.lk.T,  mis  am 
could  then,  i  n  the  absence  of  mixing,  r  ,  s  e  to  t  r.e  1 5  0  ms  .eve.. 
Winds  throughout  the  troposphere  were  very  light  from:  oemeer.  Ilf 
and  2 "  0 0  and  at  no  height  exceeded  5  m  s-'-.  The  freezing  level 
was  at  3.3km,  with  unmixed  cl  oudy  air  perhaps  1°  warmer.  :  r.  t  r.  e 
afternoon  a  few  vigorous  showers  developed. 

Two  weak  isolated  echoes  30  and  40km  East  of  Chiiboltor.  are 
clearly  visible  in  Figure  4  which  is  a  2°  e leva- i on  PPI.  A 
further  ten  kilometers  distant  at  50  k  m  range  is  a  much  more 
extensive  row  of  vigorous  cells  having  reflectivities  above 
50dBZ,  This  PPI  shows  that  at  an  altitude  of  1.4km  an  isolated 
echo  at  40  km  range  and  5km  lateral  distance  has  areas  where  2  is 
20  to  3 OdB 2  accompanied  by  anomalously  high  2DP.  values  aoove  3dB. 
This  combination  is  so  far  from  the  predictions  of  a  Marshall- 
Palmer  raindrop  size  distribution  (see  Figure  3)  that  42  PHI  and 
3  PPI  scans  were  made  through  this  ceil  in  30  minutes  in  o  r  a  e  -  to 
trace  its  evolution;  each  RHI  being  separated  by  only  1°,  wmch 
is  equivalent  to  a  distance  of  700m  at  this  range.  The  isolated 
cell  at  30km  range  and  -6km  lateral  distance  has  2  values  above 
20dBZ  but  is  not  accompanied  by  any  unexpectedly  high  values  of 
z  D  R  ■ 

The  first  vertical  section  through  the  most  intense  part  of 
this  weak  echo  obtained  at  1635-40,  three  minutes  before  the 
PPI,  is  displayed  in  Figure  5.  A  vertical  profile  of  Z  and  ZDR 
at  the  gate  where  the  maximum  Z  occurs  is  plotted  in  Figure  6. 
RHIs  at  neighbouring  azimuths  show  lower  2  values.  These  data 
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confirm  that  the  cell  had  a  Z  which  only  locally  reached  3  4  d  B  " 
whereas  the  value  of  ZDR  widely  exceeds  2dB  with  a  maximum  of 
4 . 1  d B  where  Z  is  less  than  30dBZ.  In  Figure  5  the  high  Z  q  R 
region  extends  in  a  narrow  column  almost  to  the  ground.  Half- 
melted  ice  particles  might  give  such  a  large  ZDR  but  they  would 
need  to  be  large  to  survive  such  a  fall  without  complete  melting, 
and  raindrops  would  need  to  have  a  diameter  of  about  5  or  6mm  to 
assume  such  an  oblate  shape.  in  either  case  the  low  values  of  Z 
must  imply  that  such  large  particles  are  present  in  very  low 
concentrations  . 

This  cell  subsequently  intensified  rapidly,  but  because  of  the 
light  westerly  winds  did  not  change  azimuth  appreciably  and  the 
range  only  increased  by  about  1km  every  10  minutes  or  so.  The 
actual  3  OdB  Z  echo  in  Figure  5  did  not  grow,  but,  perhaps  not 
suprisingly  in  view  of  the  large  particles  below  the  maximum 
echo,  it  fell  to  the  ground  in  the  following  five  minutes,  as 
did  the  subsidiary  lOdBZ  echo  at  37km  range.  Instead  a  local 
maximum  in  Z  of  24dBZ  at  2.5km  altitude  acccompanied  by  a  Zqr  of 
over  3dB  located  1.4km  to  the  south  (2°  in  azimuth)  increased 
rapidly  in  intensity  and  areal  extent,  so  that  by  1641-10  the 
vertical  section  in  Figure  7  was  observed  with  the  development  of 
a  narrow  vertical  column  of  positive  ZpR.  The  vertical  profile 
displayed  in  Figure  8  taken  through  the  most  intense  part  of  this 
column  confirms  values  of  Zqr  of  4dB  above  the  freezing  level 
with  Z  reaching  38dB  at  4.8km  altitude.  Two  minutes  later  the 
vertical  sections  indicate  little  change  in  the  column  but  show 
that  the  maximum  height  of  the  20dBZ  contour  has  risen  by  just 
over  1km,  equivalent  to  a  vertical  velocif'  of  almost  10m  s--1 
A  typical  RHI  obtained  a  further  five  minutes  later  at  1  648-  10  is 
displayed  in  Figure  9  and  shows  a  quite  different  character. 
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Although  the  cloud  has  grown  further  and  at  3km  height  .1  r.  as 
reached  5  4  d  B  Z  ,  a  traditional  hail  signature,  the  column  of 
positive  zDR  has  vanished.  The  vertical  profile  of  Z  and  ZDfj 
plotted  in  Figure  10  highlights  this  transformation.  the  Z  p  p 
signature  is  now  typical  of  that  commmoniy  observed  in  mature 
glaciated  cumulonimbus  clouds.  Above  the  freezing  level  Znp  is 
generally  within  O.ldB  of  zero,  this  is  usually  interpreted  as 
tumbling  hail  or  graupel  pellets.  Below  the  zero  degree  isotherm 
2 q p  gradually  increases  frc-  0  to  5dB  over  a  distance  of  2km, 
consistent  with  the  particles  slowly  melting  to  form  large  oblate 
raindrops  near  to  the  ground.  Three  minutes  later  Z  values  of 
61dB2  were  observed.  The  Z[>p  gradient  in  Figure  10  snouid  be 
contrasted  with  the  very  abrupt  fail  off  in  3  R  above  5km  :  r, 
Figure  8,  which  is  difficult  to  explain  in  terms  of  partic.es 
falling  and  melting.  A  more  plausible  interpretation  is  to 
consider  the  sharp  change  occurring  at  the  boundary  where  a 
rising  parcel  of  air  containing  supercooled  raindrops  penetrates 
upwards  into  the  glaciated  region  of  the  cloud. 

The  evolution  from  tne  first  RHI  in  Figure  5  to  the  ma-ure 
cloud  in  Figure  9  took  only  12  minutes  with  the  e c no  top  rising 
nearly  4km  at  an  average  velocity  of  5  m  s'*.  The  two  minute 
period  when  there  was  a  column  of  positive  Z^p  extending  towards 
the  - 1 0 0 C  isotherm  appeared  to  coincide  with  the  most  vigorous 
growth  with  the  Z  echo  rising  at  about  10m  s "  -  .  It  is 
interesting  to  note  that  the  second  isolated  weak  Z  echo  m  tne 
PPI  of  Figure  4  which  is  at  a  range  of  30km  and  a  lateral 
distance  of  -6km,  which  was  not  accompanied  by  any  large  positive 
ZDp,  failed  to  develop. 

3(b)  6  July  1606  to  1615 

The  1200  sonde  ascent  from  Crawley  indicated  a  similar 
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structure  to  that  already  described  on  the  7  July.  The 
calculated  freezing  level  was  3.2km  and  the  winds  again  were  less 
chat  5  m  s-3,  but  from  a  direction  between  170  and  180°.  An 
isolated  daughter  cell  was  first  detected  in  the  5°  elevation  PPI 
which  is  displayed  in  Figure  11.  This  cell  is  at  30km  range  and 
7km  lateral  distance,  at  larger  lateral  distances  no  returns 
were  detected,  but  at  2km  lateral  distance  only  3km  from  the  cell 
is  the  start  of  a  long  narrow  band  of  intense  echoes  extending 
aout  50km  to  the  south.  No  echo  was  detected  at  the  position  of 
the  isolated  cell  in  the  4°  elevation  PPI  where  the  height  of  the 
beam  at  this  range  would  be  about  2km.  Because  the  5°  elevation 
showed  such  anomalously  high  ZDR  values  of  2.4dB  accompanied  by  a 
Z  echo  of  2  5dBZ ,  a  series  of  25  RHI  scans  through  the  cell  was 
made  in  the  12  minutes  following  the  PPI.  Each  RHI  was 
separated  by  2°,  equivalent  at  this  range  to  a  horizontal 
distance  of  1km.  As  the  southerly  wind  was  so  light,  the  echo 
was  remained  between  azimuths  141  and  149  throughout  this  period, 
advecting  by  less  than  1km  in  20  minutes. 

Fig  12  displays  the  contours  in  Z  and  ZpR  f°r  the  first  RHI 
at  1607-00  through  the  most  intense  section  of  this  daughter 
cell.  The  daughter  cell  has  a  maximum  Z  of  33dBZ  but,  as  in  the 
example  of  7  July,  is  accompanied  by  much  nigher  values  of  ZDR 
than  would  be  predicted  by  most  commonly  observed  raindrop  size 
distributions.  This  is  further  demonstrated  by  the  vertical 
profiles  in  Z  and  ZDR  in  Figure  13,  where  values  of  ZDR  of  3dB 
are  found  where  Z  is  only  30dBZ.  By  similar  arguments  to  those 
employed  in  Section  3(a)  these  particles  must  be  large  raindrops 
or  large  melting  ice  particles  and  so  the  low  Z  value  implies  the 
particles  are  present  in  very  low  concentrations. 

This  daughter  cell  grew  rapidly  and  the  vertical  section  in 
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Figure  14  was  obtained  6  minutes  later  at  1613-20  with  the 
highest  2  value  of  51dBZ  coinciding  with  the  maximum  observed  ZDR 
of  5.1dB.  Again  the  Z^p  contours  show  a  well  developed  narrow 
vertical  column  containing  the  highest  values  of  differential 
reflectivity  which  extends  over  1.5km  above  the  freezing  level  to 
where  the  temperature  is  about  -10°C.  As  this  feature  is 
virtually  absent  on  the  neighbouring  RHIs  we  conclude  that  the 
diameter  of  the  column  is  only  about  1km.  An  RHI  2  minutes 
later  still  shows  the  Zqr  column,  while  Z  has  further  increase  to 
61dBZ  and  the  echo  top  has  risen  a  further  1  km  implying  an 
upward  velocity  of  close  to  10  m  s'l. 

The  evolution  is  simpler  than  the  previous  example  in  that 
only  one  maximum  in  the  echo  could  be  distinguished,  however  by 
the  time  of  Figure  14  the  anvil  from  the  band  of  more  active 
cells  visible  in  the  earlier  PPI  is  only  2km  away  from  the 
daughter  ceil.  The  maximum  echo  in  this  anvil  is  about  20dBZ, 
and  Figure  15  shows  that  by  1618-20  it  had  merged  with  this 
developing  cell  which  now  has  a  maximum  echo  of  65dBZ.  The  ZpR 
in  this  Figure  is  quite  revealing.  At  3  0 km  range  the  positive 
ZDR  is  confined  to  regions  below  the  freezing  level  indicating 
ice  above,  but  at  33km  the  positive  ZDR  extends  to  a  height  of 
4km  with  a  maximum  of  5dB  at  2.3km  altitude;  above  this  there  is 
a  region  of  negative  Zqr  with  a  peak  value  of  -1.4dB  occurring  at 
5km  where  Z  is  62dBZ  .  These  various  regimes  of  ZDR  can  be 
clearly  identified  in  the  vertical  profile  of  Z  and  ZDR  passing 
through  this  region  of  negative  Zdr  which  is  displayed  in  figure 
16.  Above  6km  Z  is  high  and  ZDR  is  very  close  to  zero,  and  this 
probably  results  from  tumbling  graupel  and  hail  pellets,  below 
the  freezing  level  (3.2km!  ZDR  increases  rapidly  and  this  is 
consistent  with  large  raindrops  or  large  oblate  melting 
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hydrometeors . 

Negative  values  of  ZqR  such  as  those  between  3.2  and  6km  in 
figure  16  have  not  been  previously  observed  with  the  Chilbolton 
radar.  In  Section  4  we  will  analyse  all  possible  effects  which 
could  give  rise  to  artificial  spurious  values  of  negative  z  p  R , 
but  if  for  the  moment  we  assume  that  the  data  is  valid,  then  the 
ony  plausible  explanation  is  that  this  negative  value  is  due  to 
oblate  or  concical  hailstones  falling  with  some  degree  of 
alignment  so  that  on  the  average  the  vertical  diameter  is  about 
20%  greater  than  the  vertical.  This  feature  is  very  localised 
and  quite  transitory.  No  negative  values  of  ZDR  can  be  found  at 
the  neighbouring  azimuths  only  1km  away,  and  in  the  previous  RHI 
at  the  same  azimuth  as  Figure  16  there  is  just  one  gate  showing  a 
negative  value  of  -0.7dB. 

The  major  evolution  appears  to  be  similar  to  the  7  July 
case  previously  described  ,  in  that  a  cell  having  low  Z  but  high 
ZqR  grew  very  rapidly.  The  highest  value  of  Z  observed  in  the 
cell  increased  from  31dB  to  65dBZ  in  only  12  minutes. 

3(c)  22  August  13.58  to  14.14 

The  evolution  of  an  isolated  shower  on  this  day  is  shown  in 
Figures  17,18  and  19.  Sonde  ascents  showed  that  the  wind  was 
again  very  light  at  all  levels  and  that  the  freezing  level  was 
about  3.1km.  The  convection  was  of  quite  limited  depth  with  the 
20dBZ  echo  contour  never  rising  above  6.6km  altitude,  but  in 
spite  of  this  the  plots  show  the  development  of  an  intense  echo 
of  over  60dBZ  which  is  associated  with  a  ZpR  of  8.3dB,  the 
highest  yet  recorded  with  the  Chilbolton  radar. 

The  first  RHI  through  the  most  intense  part  of  the  isolated 
developing  cell  obtained  at  1358-40  is  displayed  in  Figure  17  and 
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shows  a  large  area  above  the  freezing  level  where  values  of  Z3P 
exceed  4dB  and  Z  is  40dBZ.  The  echo  top  (20dBZ)  is  at  5.6km,  but 
above  5km  values  of  ZqR  are  essentially  zero.  Although  the  high 
values  of  Z^p  extend  down  to  the  ground  the  values  of  Z  are  much 
less  at  these  lower  altitudes,  indicating  that  the  major  part  of 
the  precipitation  is  still  above  the  zero  degree  isotherm  but 
that  a  few  large  precipitation  particles  are  reaching  the  ground. 
Six  minutes  later,  as  shown  in  Figure  18,  the  values  of  Z  in  the 
centre  of  the  cell  reach  over  63dBZ  but  the  ecno  top  at  this 
azimuth  is  only  5km.  These  high  values  of  Z  were  accompanied  by 
a  ZpR  value  of  8.3dB  at  one  gate.  This  gate  was  at  2.8km 
altitude,  very  close  to  the  freezing  level,  and  was  situated 
in  the  centre  of  a  region  900m  wide  and  600m  deep  where  Zpp  was 
above  6dB. 

On  this  day  the  winds  were  very  light  and  ,  in  the  absence 
of  shear,  the  intense  echo  in  Figure  18  fell  to  ground  and  the 
cell  decayed.  The  compact  region  of  ZDP  above  6dB  fell  1km  by  t.ne 
time  of  the  next  RHI  140  seconds  later,  and  after  a  further  140 
seconds  it  was  just  discernible  1.2km  lower  at  an  altitude  of 
only  600m.  By  this  time  ZDR  was  everywhere  zero  above  the 
freezing  level,  with  positive  values  confined  to  altitudes  below 
3km.  This  is  the  pattern  usually  observed  with  mature  convective 
storms.  Figure  19  was  obtained  eleven  minutes  after  Figure  18, 
and  shows  that  the  highest  values  of  ZDR  were  only  2dB  where  Z 
was  above  50dBZ  implying  raindrop  size  distributions  much  closer 
to  the  Marshall-Palmer.  From  Figure  3  if  Nc  is  8000  m' 3  mm'6 
then  the  contours  of  ldB  and  2dB  in  ZqR  should  be  near  to  the 
contours  of  35dBZ  and  48dBZ  respectively.  This  approximates  to 
the  situation  in  Figure  11. 
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4.  POSSIBLE  ARTIFACTS  IN  DIFFERENTIAL  REFLECTIVITY  Y  EASY  PE. v  ENTS. 

Before  embarking  upon  an  interpretation  of  the  rather 
unexpectedly  nigh  values  of  LDR  reported  in  r he  aoove  convective 
clouds,  we  consider  whether  these  ODservat  lor.s  could  have  been 
influenced  by  any  possible  spurious  factors.  Beard  et  al  (198?) 
have  suggested  that  raindrops  oscillate  and  that  their  mean  shape 
is  slightly  more  spherical  than  the  equilibrium  one;  if  tnis 
effect  was  significant  then  a  given  value  of  2DR  should 
correspond  to  an  even  larger  raindrops  size  than  that  assumed  in 
this  paper. 

We  now  examine  the  effect  of  side  lobes  in  the  presence  of 
reflectivity  gradients.  The  symbols  superposed  on  the  five 
Figures  showing  the  vertical  profiles  represent  each  separate 
data  point  recorded  for  every  quarter  degree  increment  m  antenna 
elevation.  This  angular  resolution  is  equal  to  the  beam width 
of  the  Chilbolton  antenna.  In  all  these  profiles  changes  in 
reflectivity  of  20dBZ  never  occur  over  such  small  angular 
changes,  so  for  the  Chilbolton  antenna,  ihe  signal  detected  by 
the  side  lobes  is  always  much  less  than  that  sensed  by  the  main 
beam.  Consequently  the  contribution  to  the  measured  differential 
reflectivity  by  any  mismatched  polarisation  characteristics  of 
the  sidelobes  is  negligible.  It  is  also  clear  from  Figure  16 
that  the  negative  values  of  ZDR  coincide  with  the  maximum  values 
of  Z  of  over  60dBZ.  Similarly  the  maximum  values  of  Z  D  R  in 
Figure  12  occur  where  Z  attains  its  highest  value.  If  these 
features  in  ZDR  were  artifacts  due  to  side  lobes  then  they  should 
be  found  where  Z  is  changing  rapidly. 

Differential  attenuation  can  result  in  an  apparent  reduction 
in  the  values  of  ZDR-  If  large  oblate  precipitation  particles 
lie  between  the  cell  and  the  radar  antenna,  then  the  incident 
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horizontal  radiation  will  be  attenuated  by  more  than  the  vertical 
and,  consequently,  the  received  power  reflected  even  from 
spherical  particles  will  be  less  in  the  horizontal  than  in  the 
vertical,  leading  to  an  apparent  negative  value  of  Zqr.  In 
Figure  14,  where  we  have  tentatively  ascribed  the  negative  ZDR 
to  conical  hail,  there  is  no  echo  bet  ween  the  cell  and  the 
antenna  and  so  differential  attenuation  cannot  be  responsible. 
It  is  worth  commenting  that  even  using  10cm  radiation  this 
effect  can  occasionally  be  important,but  is  a  sensitive  function 
of  the  oblate ness  of  the  intervening  precipitation.  Calculations 
assuming  an  exponential  raindrop  size  distribution  show  that  for 
a  rainfall  of  10  0mm /hr  having  a  Z D R  of  ids  the  differential 
attenuation  is  only  O.OldB/km,  but  if  2  DR  is  3dB  t  hen  the 
differential  attenuation  is  increased  tenfold.  These  figures 
are  an  order  of  magnitude  higher  at  5.6cm.  Negative  values  of 
Z DR  of  below  -0.2dB  can  be  detected  in  the  PPI  in  figure  11,  at  a 
range  of  29  to  30km  and  a  lateral  range  of  -7km;  the  lowest 
value  recorded  being  -0.6dB  at  one  gate.  The  radiation  between 
this  region  and  the  antenna  has  traversed  a  2km  region  having  a  2 
of  over  50  d B 7.  and  a  differential  reflectivity  above  3dB.  A 
differential  attenuation  of  O.ldB- km  can  account  for  the  slightly 
negative  values  of  ZDR.  In  none  of  the  other  sections  discussed 
is  this  effect  important. 

The  observation  in  mature  intense  convective  storms  that  the 
values  of  Zqr  are  so  close  to  zero  for  all  values  of  Z  up  to 
65dBZ  gives  us  additional  confidence  that  there  is  no  systematic 
bias  between  the  processing  of  the  signals  of  the  vertical  and 
horizontal  polarisations.  We  conclude  that  the  reported 
differential  reflectivity  observations  are  valid. 
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5. 


INTERPRETATION  OF  THE  1933  RESULTS. 


The  three  evolutions  described  in  Section  3  are  the  examples 
where  the  echoes  were  located  and  observed  at  an  early  stage  :n 
their  growth.  The  most  striking  aspects  of  the  evolution  are  the 
occurrence  of  high  positive  values  of  Z^r  accompanied  by  fairly 
modest  values  of  z  as  shown  in  Figures  5,  12  and  17.  in  view  of 
the  arguments  advanced  in  Section  2  we  believe  that  these  high 
values  of  Zqr  imply  large  particles.  In  the  early  stages  the 
high  values  of  ZDR  occur  both  below  and  above  the  zero  degree 
isotherm  indicating  that  they  are  probably  large  raindrops.  If 
they  were  due  to  wet  ice  particles  either  in  wet  growth  or 
melting,  then  they  would  not  be  expected  to  be  more  oblate  than 
a  raindrop  of  the  corresponding  size.  We  shall  offer  an 
explanation  assuming  that  these  high  values  of  Z  r>  R  are  due  to 
large  raindrops,  but  shall  retain  the  possibility  that  they  are 
due  to  even  larger  wet  ice  particles. 

Figure  5  and  the  accompanying  vertical  profile  in  Figure  6 
show  that  in  this  first  RHI  of  the  7  July  example,  Z  values 
near  to  30dBZ  were  accompanied  by  ZDR  of  about  3dB.  Assuming  a 
monodispersed  distribution  of  raindrops  then  Figure  2  indicates 
that  the  cloud  contains  drops  of  about  4mm  diameter  present  in  a 
concentration  of  rather  less  than  one  per  cubic  meter.  The  more 
realistic  exponential  distribution  plotted  in  Figure  3  using 
Pruppacher  and  Pitter  drop  shapes  predicts  a  value  of  N0  of  about 
8  (instead  of  8000  in  the  normal  Marshall  Palmer  distribution) 
and  a  rainfall  rate  of  0.2  mm  'hr.  Even  with  the  more  oblate  drop 
shapes  of  Cooper  et  al  N0  is  40  and  R  is  0.35  mm/hr.  The  first 
RHI  obtained  in  the  6  July  example  in  Figure  12  and  the  profile 
in  Figure  13  again  show  a  Z  of  slighlty  less  than  30dBZ  where  ZDr 
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is  ]ust  above  3dB.  These  observations  again  suggest,  using  either 
the  monodispersed  or  exponential  model  that  there  is  a  very  low 
concentration  of  large  (above  4mm)  supercooled  raindrops  present 
in  the  cloud.  In  this  case  the  rainfall  rate  is  considerably 
less  than  1mm /hr,  whereas  an  empirical  Z-R  relationship  would 
indicate  about  2, 5mm, 'hr..  The  example  on  the  22  August  shows  a 
similar  deviation  from  the  Marshall-Palmer  distribution.  In  the 
centre  of  the  echo  is  a  large  area  having  Z  values  of  about  40dBZ 
occurring  where  ZDR  is  about  4dB,  suggesting  a  mean  raindrop  size 
of  5  or  6mm.  Figure  3  shows  that  this  data  cannot  be  fitted  to 
an  exponential  raindrop  size  distribution  if  the  Pruppacher  and 
Pitter  drop  shapes  are  used,  but  a  fit  is  possible  using  the  more 
oblate  shapes  of  Cooper  et  al,  and  implied  values  of  N0  are  at 
least  two  orders  of  magnitude  lower  than  the  Marshall-Palmer 
value . 

The  sizes  and  concentration  of  raindrops  calculated  above 
relate  to  the  region  of  the  cloud  where  the  reflectivity  is  a 
maximum;  in  other  parts  of  the  cloud  the  implied  concentrations 
are  even  lower.  Evaporating  ra  inshafts  could  lead  to  small 
concentrations  of  large  drops,  but  these  should  be  fairly  local 
phenomena  not  occurring  near  the  reflectivity  maximum  of  a 
developing  cloud.  Such  a  low  concentration  of  large  drops  would 
be  difficult  to  confirm  with  aircraft  because  the  cloud  is 
evolving  very  rapidly  and  conventional  probes  would  not  obtain  a 
statistically  significant  sample  in  a  single  1  or  2km 
penet  ration. 

The  evolution  which  emerges  from  this  limited  study  is  rather 
unusual.  A  small  number  of  large  (4mm)  drops  form,  which  are 
liquid  even  at  -10°C.  Rapid  intensification  of  the  echo  follows, 
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with  Z  rising  at  a  rate  of  up  to  30dBZ  in  ten  minutes,  with 
subsequent  glaciation  and  disappearance  of  positive  values  of  Zgp 
aoove  the  zero  degree  isotherm.  The  most  rapid  upward  growth  of 
the  echo  at  a  velocity  of  about  10m  s"1  appears  to  coincide  with 
partial  glaciation,  although  values  of  ZDp  are  mostly  zero  where 
the  temperature  is  below  freezing,  there  is  a  small  column  of 
positive  Zopj  extending  to  an  altitude  where  the  temperature  is 
about  -10°C.  This  column  is  a  transitory  phenomenon  persisting 
for  less  than  ten  minutes;  it  is  of  circular  cross  section  with  a 
diameter  of  between  1  and  2km,  but  the  values  of  typically 
reach  4dB  within  the  column.  As  the  cell  1  r  -  s  its  vigour  the 
echo  descends  and  the  area  of  high  a. so  subsides  towards  the 
ground  without  any  very  marked  change  in  its  magnitude. 

It  seems  reasonaDle  that  a  rapidly  rising  parcel  of  moist 
air  could  rise  to  an  altitude  where  the  temperature  is  -10°C 
before  any  of  the  surrounding  ice  particles  mix  in  and  cause 
rapid  glaciation.  It  could  be  argued  that  large  ice  particles  in 
wet  growth  can  account  for  these  observations,  although  their 
presence  would  be  equally  difficult  to  explain  at  these  early 
stages  of  cell  growth,  but  as  suggested  before,  the  surplus 
water  should  be  shed  and  not  lead  to  large  positive  2  D  r  . 
Furthermore  as  the  particles  fall  and  melt,  the  value  of  Zqr 
should  rise  even  further  and  this  is  not  observed. 

Once  the  clouds  in  this  study  glaciate  the  ice  particle 
spectrum  will  intially  be  similar  to  the  raindrop  spectrum. 
Federer  and  Waldvogel  (1975)  report  mean  spectra  for  hail  of  size 
less  than  22mm  having  an  average  N0  of  12.1m“3  m'1  and  a  D0  of 
8.73mm.  Similar  values  of  D0  and  N0  are  apparent  from  Figure  3 
for  the  initial  raindrop  spectra  in  these  three  examples.  There 
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is  some  evidence  from  aircraft  penetrations  of  developing  cells 
(Musil  et  al,  19761  that  the  value  of  N0  may  be  as  low  as  10  or 
100  m"3.  : ;  even  smaller  concentrations  of  large  supercooled 
raindrops  do  exist  in  the  early  stages  of  precipitation 
development  then  they  would  be  the  natural  favoured  embryos  which 
could  then  grow  rapidly  into  hail  without  exhausting  the  supply 
of  liquid  water.  Consequently  the  requirement  of  embryo  sorting 
(Browning,  1978)  for  hail  growth  may  not  be  quite  so  stringent, 
"acd.t  et  ai  (I960)  found  that  most  hailstones  collected  in 
fn g.and  had  clear  growth  centres  suggesting  large  frozen  drop 
emor yos . 

In  one  example  we  infer  that  large  hail  causes  a  negative 
va.ue  of  although  this  seems  to  be  a  rare  occurrence.  More 
normally  the  value  of  ZDR  is  nearer  to  zero,  indicating  that  hail 
‘umbles  as  it  falls  or  is  asymmetric.  A  much  more  reliable  hail 
signature  at  low  altitudes  in  the  cloud  is  a  region  of  zero  Z^p 
where  Z  is  very  high,  instead  of  the  more  normal  postive  in 
neavy  rain  (Illingworth  et  al,  1986). 

6.  CONCLUSION. 

From  a  study  of  a  number  of  vigorous  isolated  convective 
storms  over  the  UK  during  the  summer  of  1983  we  draw  the 
following  tentative  conclusions. 

a)  Isolated  weak  echoes  having  moderate  Z  but  accompanied  by 
large  positive  differential  reflectivity  (  Z  D  R  1  are  1  1  *  e  1  ^  t0 
intensify  very  rapidly. 

b)  These  anomalous  low  z/  high  positive  ZqR  first  echoes 
indicate  very  low  concentrations  of  large  precipitation 
particles.  Although  we  cannot  rule  out  the  possibility  that  they 


3  re  1  a  r  9  e  ice  p  a  r  1 1  c  1  e  s  in  wet  growth  this  s  e  e  m  s  most  unlikely, 
and  we  propose  than  they  are  large  ( "4mm/  supercooled  raindrops. 
Present  t.ieeries  of  precipitation  development  do  not  predict  the 
appearance  of  tnese  large  drops  at  such  an  early  stage. 

c)  Such  low  concentrations  of  large  drops  could  be  very 
efficient  nail  embryos  growing  to  large  hailstones  after  freezing 
because  of  the  lack  of  competiition  for  the  cloud  water.  This 
may  relax  requirements  for  embryo  sorting. 

d)  Narrow  columns  of  positive  ZDR  extending  by  up  to  2km  aoove 
the  zero  degree  isotherm  occur  during  vigorous  convection,  when 
the  echo  top  is  rising  rapidly  and  tend  to  preceed  high,  values  of 
Z .  They  are  generally  circular  m  cross  section  of  diameter 
between  1  and  2km,  the  remainder  of  the  cloud  at  these  heights 
having  ZDp  values  very  close  to  zero.  We  believe  that  these 
columns  may  be  due  to  large  supercooled  raindrops  in  an  ascending 
updraught,  and  that  their  subsequent  disappearance  indicates 
rapid  glaciation. 

Further  data  sets  obtained  with  me  Chilbolton  radar 
should  provide  more  evidence  on  the  size,  concentration  and  types 
of  hydrometeors  during  the  evolution  of  convective  clouds.  The 
measurement  of  additional  polarisation  parameters  should  provide 
definite  confirmation  of  the  phase  of  the  scattering  particles, 
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LEGENDS 


Figure  1:  The  values  of  Zqr  predicted  from  Gans  scattering  theory 
for  oblate  spheroids  as  a  function  of  their  minor  to  major  axis 
ratio,  (a/b).  Figure  axis  vertical.  Water  -  curve  1 ;  solid  ice 
of  0.92  g  cm~3  -  curve  2;  ice  air  mixture  of  0.5  g  c  m  ~  3  .  curve 
3;  and  an  ice  air  mixture  of  O.lg  cm  -  3  -  curve  4. 

Figure  2:  The  variation  of  reflectivity,  Z,  with  differential 
refectivity,  Zqr,  for  a  monodispersed  distribution  of  raindrops 
having  equivalent  diameter  D.  R:  Rainfall  rate  of  1  mm  hr'l. 

N :  drop  concentrations  of  100  and  1000  m " 3 .  Solid  curves 
Pruppacher  and  Pitter  shapes,  dotted  curves  shapes  extrapolated 
from  Cooper  et  al. 

Figure  3:  The  variation  of  Z  with  Z  p  r  predicted  for  the 
exponential  distribution  of  raindrops  in  Equation  <21.  Solid 
curves  for  N0  =  8000  m~3  mm"'-  and  R  =  1  mm  hr~^  using  Pruppacher 
and  Pitter  shapes,  dotted  branches  for  the  more  oblate  shapes  of 
Cooper  e t  a  1 . 

Figure  4:  A  horizontal  section  (ppi)  at  2°  elevation  at  1638-00 
on  7/7/83  showing  two  isolated  cells.  In  this  and  the  following 
Figures  the  solid  contour  lines  in  Z  are  drawn  every  lOdBZ 
starting  at  lOdBZ,  and  in  Z^r  every  ldB  starting  at  ldB.  The 
broken  contour  line  in  both  the  Z  and  the  Zqr  plots  is  for  Z=5dBZ 
and  serves  to  delineate  the  extent  of  the  cloud. 

Figure  5:  A  vertical  section  at  1635-40,  azimuth  93°,  through  the 
most  intense  part  of  the  echo  at  40km  range  in  figure  4.  Contours 
of  Z  and  Zqr  as  in  Figure  4 

Figure  6.  A  vertical  profile  of  Z  and  ZDB  at  40.8km  range,  the 
gate  where  Z  is  a  maximum  in  Figure  5.  Below  3km  altitude  values 
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of  ZqR  are  anomalously  high. 

Figure  7.  An  RHI  at  1641-10  through  the  cell  in  Figure  5  at 
azimuth  95°  showing  the  development  of  a  narrow  column  of 
positive  ZqR  above  4dB  which  extends  almost  2km  above  the 
freezing  level. 

Figure  8:  A  vertical  profile  of  Z  and  Z^p,  at  range  39.9km  through 
the  column  of  positive  Z[)R  in  Figure  7. 

Figure  9:  An  RHI  obtained  at  1648-10  and  azimuth  94°  through  the 
same  cell  as  the  previous  five  figures,  but  now  Z  and  Zqr  are 
more  typical  of  a  mature  glaciated  cloud,  with  ZDR  only  positive 
below  the  freezing  level. 

Figure  10:  A  vertical  profile  of  Z  and  ZqR  at  a  range  of  40.8km 
in  Figure  9.  Z  p  R  values  are  zero  above  3km  consistent  with 
tumbling  graupel,  but  become  increasingly  positive  below  3km  as 
the  particles  melt. 

Figure  11:  A  PPI  obtained  at  1606-00  on  6/7/83  showing  an 
isolated  cell  with  unexpectedly  high  ZDR  at  30km  range  and  7km 
lateral  distance,  separated  by  3km  from  an  extensive  area  of 
intense  echoes.  An  area  of  negative  ZDR  at  29km  range  and  lateral 
distance  -7km  (marked  with  an  N  )  is  believed  due  to  differential 
attenuation  and  is  discussed  >n  Section  4. 

Figure  12:  An  RHI  at  1607-00,  azimuth  148°,  through  the  most 
intense  part  of  the  isolated  cell  in  figure  11.  High  positive 
ZqR  accompanies  a  weak  Z. 

Figure  13:  A  vertical  profile  of  Z  and  ZDR  at  32.1km  range 
through  the  maximum  z  echo  in  Figure  12.  Comparison  with  Figure  3 
confirms  the  anomalously  high  values  of  ZDR. 

Figure  14:  An  RHI  at  1613-20,  azimuth  146°,  showing  the 
development  of  a  narrow  column  with  ZDR  values  exceeding  5dB 
which  extends  over  1.5km  above  the  freezing  level. 
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Figure  15:  An  RHI  at  1618-20,  azimuth  145°,  through  the  same 
cell  as  the  previous  four  Figures.  An  additional  -0.5dB  contour 
in  Zjjr  is  plotted  (marked  with  an  N),  and  occurs  where  Z  is  above 
60dBZ.  The  negat i ve  Zqr  is  believed  to  be  due  to  conical  hail. 
Figure  16:  A  vertical  profile  of  Z  and  Zqr  at  33km  range  through 
the  negative  Zqr  region  in  Figure  14,  demonstrating  that  the 
negative  Zqr  occurs  where  Z  is  a  maximum. 

Figure  17:  The  first  RHI  through  a  developing  cell  at  1358-40  on 
22/8/83,  azimuth  78°.  Note  the  high  Z^r  above  the  freezing 
level  . 

Figure  18:  An  RHI  through  the  cell  in  figure  17  obtained  at 

1404.40,  azimuth  78°.  The  highest  Zqr  contour  is  7dB  with  a 
maximum  value  at  one  gate  of  8.3dB.  Echo  top  is  only  5km. 

Figure  19:  An  RHI  obtained  11  minutes  later  than  Figure  18  at 

1415.40,  azimuth  80°.  The  large  precipitation  particles  have 
fallen  to  ground  and  Z  and  Zqr  are  much  closer  to  those  expected 
for  a  Marshall-Palmer  drop  size  distribution. 
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